Conduction mechanisms in biphenyl-dithiol single-molecule junctions 
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Based on density-functional theory calculations, we report a detailed study of the single-molecule 
charge-transport properties for a series of recently synthesized biphenyl-dithiol molecules [D. Von- 
lanthen et al, Angew. Chem., Int. Ed. 48, 8886 (2009); A. Mishchenko et al, Nano Lett. 10, 156 
(2010)]. The torsion angle ip between the two phenyl rings, and hence the degree of tt conjugation, 
is controlled by alkyl chains and methyl side groups. We consider three different coordination ge- 
ometries, namely top-top, bridge-bridge, and hollow-hollow with the terminal sulfur atoms bound 
to one, two, and three gold surface atoms, respectively. Our calculations show that different coor- 
dination geometries give rise to conductances which vary by one order of magnitude for the same 
molecule. Irrespective of the coordination geometries, the charge transport calculations predict a 
cos^ Lp dependence of the conductance, which is confirmed by our experimental measurements. We 
observe that the calculated transmission through biphenyl dithiols is typically dominated by a single 
transmission eigenchannel formed from tt electrons. Only for a single molecule with a completely 
broken conjugation we find a perfect channel degeneracy for the hollow-hollow-type contact in our 
theory. 

PACS numbers: 85.65. +h, 73.63.Rt, 73.23.Ad, 31.15.es 



I. INTRODUCTION 

After the first realizations of single- or few-molecule 
contacts,^^ a major theme of research represents the 
controlled fabrication of molecular junctions with desired 
properties. As an example, several groups demonstrated 
recently that the conductance of junctions containing 
biphenyl derivatives can be controlled by the torsion an- 
gle tp between the two phenyl rings.*"* However, not 
only the investigation of the conductance of more com- 
plex molecules has become feasible j^^M^ but also the study 
of additional aspects such as the signature of molecu- 
lar vibr ation s in the electric currentjl^Sl current-induced 
heatingjlIiEll and the distinction of electron or hole con- 
duction by measurement of the thermopower.E^l Further- 
more, information on individual conduction channels can 
be obtained by use of superconducting electrodes^^ or 
shot-noise measurements.!^ All these advances allow for 
a better characterization of the single-molecule charge 
transport when compared with theory. 

Despite experimental and theoretical achievements, 
measurement and modeling of electron transport in 
molecular junctions are still challenging tasks. This is 
mainly due to the observed variability in junction con- 
ductances a nd the corresponding statistical nature of the 
experiments. '^E^'^ In this regard, calculations based on 
the approximate density functional theory (DFT) can be 
helpful to obtain a better understanding of the charge 
transport mechanisms involved and to interpret trends 
in the experimental data based on computed structure- 
transport relationships. In agreement with the experi- 



mental observations, they show, in particular, that the 
electric condu ction strongly depends on the molecular 
conformatiorPHsElHlll and the precise geometry in the 
single- molecule junctions .'^'^ 

In our recent studies we have explored the con- 
duction properties of biphenyl-dithiol (BPDT) molecules 
bound to An electrodes. For these molecules, named here 
M0-M7 and displayed in Fig. [T| the molecular conjuga- 
tion is gradually varied by the use of alkyl chains and 
methyl side groups. In this follow-up paper, we present a 
more detailed theoretical analysis of their transport prop- 
erties based on DFT calculations. We study an extended, 
systematic set of contact geometries and place special em- 
phasis on transport for perpendicular ring orientations. 
The conduction mechanisms are revealed by means of 




Figure 1: Chemical structure of the investigated molecules. 
"S" is the sulfur atom, and "R" represents the acetyl group for 
the synthesized form of the molecules, a hydrogen atom after 
the in-situ deprotection, or the Au electrode for the transport 
measurements. 
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a tight-binding model (TBM), the more frequently used 
Lorentz model (LM), and the eigenchannel decomposi- 
tion of the conductance. Our calculations suggest that 
the coordination site ("top", "bridge", or "hollow") of the 
anchoring sulfur atom at the Au surface plays a decisive 
role in conduction through the molecular junction. 

The paper is organized as follows. In Sec. [IT] we briefly 
summarize key experimental findings.^ Technical as- 
pects of the DFT and transport calculations are dis- 
cussed in Sec. Studies of the dependence of junc- 
tion conductance on the torsion angle and the substrate- 
adsorbate coordination geometry are presented and an- 
alyzed in Sec. |IV[ The main text ends with a summary 
and conclusions in Sec. |V) The appendices contain de- 
tails on the methods used to analyze the DFT results in 
Sec. |IV[ Thus, in Appendix \K\ we explain our scheme to 
construct transmission eigenchannel wavefunctions with- 
out the need to resort to Lowdin transformations when 
nonorthogonal, local basis sets are used in the DFT cal- 
culations. Finally, Appendix [B] discusses the relation be- 
tween the TBM and the LM. 



II. EXPERIMENTS 

We synthesized the BPDT molecules M0-M7 of 
Fig. [T] with acetyl-protected terminal thiol groups, 
i.e. R=C0CH3.-^^ The torsion angle ip between the 
phenyl rings is gradually varied by introduction of alkyl 
side chains of variable length or methyl groups.EI^Si This 
leaves the length of the molecule unchanged. The torsion 
angles were determined by an X-ray structure analysis of 
single crystals formed from each of these compounds, ex- 
cept for MO. 

This systematic set of molecules exhibits several re- 
markable features in single-molecule transport measure- 
ments, which we examine further below. First, stable 
junctions can be formed with gold leads by the termi- 
nal sulfur atoms after deprotection. Second, the confor- 
mation of the biphenyls is efficiently locked by the alkyl 
chains and the steric hindrance of the methyl side groups. 
Var iations of torsion angles are hence expected to be 
low Third, strongly electron-donating or electron- 
withdrawing side groups are avoided, which have been 
demonstrated to influence noticeably the single-molecule 
transport .123 

We studied the conductance of single-molecule junc- 
tions by means of a scanning-tunneling-microscopy- 
break-junction method (see Ref. 'S). The most proba- 
ble or "typical" molecular junction conductance is deter- 
mined as an average of the peak values in the conductance 
histograms, measured at three different bias voltages (see 
Fig.|2| inset). 

When we plot this typical conductance as a function of 
the torsion angle on a linear scale, we find a G = a cos^ (p 
dependence with a — 2.49 • 10~^Go, as expe cted for off- 
resonant TT-dominated charge transport .1^1^21 Here, Gq = 
2e^/h is the conductance quantum. Fig. [2] shows a semi- 
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Figure 2: (Color online) Dependence of the experimental con- 
ductance G (points) on the torsion angle ifi for M1-M7 with 
ip determined from X-ray structures. The solid line repre- 
sents a fit to G = acos'^ (p with a = 2.49 • 10~*Go. The inset 
shows conductance histograms for M7, obtained at three dif- 
ferent bias voltages. The typical single-molecule conductance 
(points in the main panel) is obtained as the average of the 
peak positions of the fitted log-normal distributions at the 
three difi'erent biases, and error bars are determined from the 
peak variations. For M7 this yields G = (9 ± 2) • IO'^Gq. 



logarithmic plot of G vs. tp. The graph reveals several 
distinct features. The conductance of Ml is lower than 
expected from the general trend. We observed a similar 
exceptional behavior of Ml in a recent investigation of 
cyano-terminated molecules (i.e. S is replaced by CN in 
Fig. l),lfil showing, however, a higher value than expected, 
and these irregularities are currently of an unclear origin. 
On the other hand, there are systematic deviations from 
the cos^ ip law for the molecules M6 and M7 with (p > 80° . 
A simple 7r-orbital modeP^El] looses its validity for large 
torsion angles, and any other than tt-tt couplings prevent 
the complete suppression of transport. We show below 
that the residual couplings are of tt-ct type. 

We observe that the conductance histograms, such 
as those shown in the inset of Fig. [2] are asymmet- 
ric with a long tail towards higher conductance val- 
ues. The broad tail region could be related to junc- 
tions with multiple molecules, modifications in substrate- 
adsorbate coordination from junction to junction, atomic 
rearrangements upon stretching, local surface roughness, 
or electrode-induced changes of the average torsion an- 
gle due charge transfer and geometric constraints. In 
contrast to results for biphenyldiamines,^ we do not ob- 
serve clear correlations between the fuU-widths-at-half- 
maximum of the conductance peaks and the expected 
differences in torsion-an gle-r elated energy barriers of the 
various BPDTs studied. SliSI jf ^ substantial part of the 
experimental conductance scatter would be due to the 
variation of the torsion angle of the biphenyl core, then 
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MO with a low energy barrier for ring rotation of about 
0.1 eV^^ should exhibit a particularly broad conductance 
distribution. However, we see no evidence to support this 
hypothesis. In agreement with conclusions from other 
works ■^^ we propose that the use of thiol anchoring groups 
leads to the variation of the single-molecule conductance 
being dominated by changes in the metal-molecule con- 
tact. 



III. THEORETICAL PROCEDURES 





R 



A. Electronic structure and geometry optimization 

Electronic structure calculations and geometry opti- 
mizations are performed within DFT. We use the quan- 
tum chemistry package TURBOMOLE 6.2.1^ For all 
calculations, we employ the standard basis set, def- 
SV(P), which is of split- valence quality with polariza- 
tion functions on all non-hydrogen atoms.l^^l^ We treat 
all molecules and contact geometries as open-shell sys- 
tems with no unpaired electro ns, an d use BP86 as the 
exchange-correlation functional.'^SllI] Total energies are 
converged to a precision of better than 10~^ a.u. and 
geometry optimizations are carried out until the change 
of the maximum norm of the Cartesian gradient is below 
10-4 a.u. 



B. Charge transport calculations 

We determine conduction properties within the 
Landauer-Biittiker formalism.**^ The energy-dependent 
transmission t(E) is expressed using standard Green's 
function techniques (see also Appendix [A| . The conduc- 
tance at low temperatures is then given by 

n 

with T„ being the transmission probability of the trans- 
mission eigenchannel n. 

In the calculations, we model the electrodes of a molec- 
ular junction as perfect semi-infinite crystals to the left 
and to the right. The molecule is connected to their sur- 
face by atomically sharp metal tips. We describe this by 
computing a finite "extended central cluster" (ECC), as 
displayed in Fig. [3] into which large parts of the metal 
electrodes are included to ensure the proper alignment of 
molecular levels with respect to Ep- 

Due to the locality of the Gaussian basis sets employed, 
we are able to partition the ECC into three subsystems, 
formed from basis states in the left (L), central (C), and 
right {R) parts. The atoms in the L and R regions of 
the ECC are assumed to represent that part of the semi- 
infinite crystal surface which couples to C. We extract 
the parameters for a description of region C and its cou- 
pling to the left and right electrode surfaces from the 



Figure 3: (Color online) Division of the ECC into the L, C, 
and 7? regions. A large number of gold atoms (around 120 in 
L and R, respectively) is used to represent the electrodes in 
the DFT calculations. 

electronic structure of the ECC. On the other hand, the 
surface Green's functions of the L or R electrodes are 
constructed using parameters obtained from a spherical 
Au cluster of several hundred atoms. This calculation 
yields a Fermi energy of Ep — —5.0 eV. With these in- 
gredients, we compute the transmission probability t(E). 
A more detailed description of our cluster-based density- 
functional approach to quantum transport can be found 
in Ref. 43, 

In transport experiments with single molecules, often 
only the low-bias conductance, proportional to the sum 
of the T„ in Eq. Q , is measured. However, also the indi- 
vidual Tn can be resolved. ^^^1^^ On the theory side, in ad- 
dition to the transmission probabilities of the conduction 
eigenchannels, also the projection of their wavefunction 
onto the central region can be obtained fr om q uantities 
at hand in the Green's-function formalism.HSHSl In order 
to construct energy-normalized transmission eigenchan- 
nel wavefunctions, which can be compared to each other, 
we proceed along the lines of Ref. |35] Our efficient pro- 
cedure, which avoids the Lowdin transformation of Ref. 
i45| is presented in detail in Appendix \K\ 

C. Contact geometries 

The statistical nature of the single-molecule conduc- 
tance experiments (see Sec. |ll]) does not provide an a 
priori assignment of representative junction geometries. 
Therefore, we have decided to study three contact struc- 
tures with different coordinations of the terminal sulfur 
atoms. The procedure adopted to determine the struc- 
ture of the ECC is summarized in Fig. [4] 

In order to model the molecular junctions, we connect 
the molecule to two Au (111) pyramids, both stemming 
from the same ideal fee Bravais lattice. We consider the 
following three types: For hollow-hollow (HH) [Fig.|4](a)] 
the S atoms of the molecule are bound at each side to 
three Au atoms, for bridge-bridge (BB) [Fig. 141(b)] to 
two Au atoms, and for top-top (TT) [Fig. |4](c)pto only 
a single Au atom. Junctions of the form HB, HT etc. 
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Figure 4: (Color online) Procedure used to set up the contact 
geometries for (a) HH, (b) BB, and (c) TT binding, respec- 
tively. 

should also occur in the experiments, but they are not 
considered here. 

For the determination of the HH and TT geometries, 
we start from the gas-phase structure of each molecule 
(with SR=H in Fig.fTl), replace the terminal H atoms by 
an S-Aui group (R=Sui in Fig.[l]), and compute ground- 
state geometries. For HH contacts, the Aui atoms are 
removed. An Auig cluster, resembling a Au (111) pyra- 
mid with a thiolated benzene attached, is computed sep- 
arately. The cluster is positioned at each side of the 
BPDT such that the S atoms on top of the pyramids co- 
incide with the S atoms of the molecule [Fig. |4ja)]. To 
obtain the TT geometries, the molecule is oriented such 
that each Aui atom coincides with a tip atom of the Au2o 
pyramids [Fig.Qc)]. To determine the equilibrium struc- 
ture for both HH and TT, the inner part is relaxed and 
only the two outermost gold layers, consisting of 6 and 
10 atoms, are kept fixed in the ideal Au fee structure. 

For the BB geometries we follow slightly different steps. 
First the terminal H atoms of the gas-phase molecules 
(with SR=H in Fig. [T]) are replaced with S, one side is 
connected to a Au2o pyramid in bridge position, while 
the other one is terminated with Aui. The outermost 
gold layers of the pyramid are kept fixed, while the rest 
is optimized. A second Au2o pyramid is finally added to 
the Aui-terminated side, where the relative distances of 
the binding S atom with respect to the new Au2o cluster 
are chosen to be the same as for the S atom in bridge 
position at the Au2o-terminated side [Fig. |4jb)]. Fixing 
again only the two outermost Au layers, the structure is 



optimized to determine the ground-state geometry. 

We note that the contact geometries do not only dif- 
fer with respect to the coordination of the sulfur atoms 
to the gold electrodes, but also in the stress exerted on 
the molecules. As visible in Fig.|4|^a), the (111) direction 
is located in the ring plane of a mono-thiolated benzene 
molecule on top of a Au pyramid. Since the S-S axis is 
along the same direction, the BPDT molecule is expected 
to adopt a minimum-energy configuration inside the HH 
junction with close to its gas-phase angle. In contrast, 
in the TT geometries the biphenyl derivative bridges the 
gold tip atoms, which are opposite to each other. In this 
case the sulfur atoms are deflected from their equilibrium 
positions, which would be located along the (111) direc- 
tion on top of the pyramids [Fig. |4|c)]. In the geometry 
optimizations we find that the phenyl ring planes of the 
biphenyl molecules tend to align parallel to the surfaces 
of the pyramids. Since this may not be possible on both 
sides of the junction, some torque is exerted. Beside ef- 
fects related to charge transfer, which may also be present 
for the HH contacts, geometric constraints thus yield an 
additional contribution to the change of the torsion an- 
gle. Similar effects as for TT are also present for the BB 
contacts, since the orientations of the phenyl rings with 
respect to the gold pyramids on both sides are generally 
different according to our construction. 

We have determined binding energies by subtracting 
the total energy of the contact geometries from those 
of the frozen separate parts, namely the left and right 
Au clusters and the S-terminated biphenyl (without hy- 
drogen on the sulfur atoms). With this procedure, we 
find the following averaged binding energies for the set 
of molecules: 5.9 ± 0.3 eV (HH), 2.9 ± 0.2 eV (BB), and 
2.2 ± 0.2 eV (TT). Hence, we find a trend of decreas- 
ing binding energies with decreasing coordination of the 
sulfur atoms to Au. 

For reasons of computational feasibility, the structural 
optimizations (and calculations of binding energies) are 
carried out with Au pyramids consisting of 19 atoms for 
HH, and 20 atoms for TT and BB, respectively. To ensure 
a proper description of the Fermi-level alignment in the 
transport calculations, the gold pyramids are extended 
to 115 (HH) and 116 atoms (BB, TT), as displayed in 
Fig. [3] All added atoms are positioned on the ideal fee 
lattice with a lattice constant a — 0.408 nm, matching 
those of the fixed layers for the smaller pyramids. No 
further geometry optimization is carried out for contacts 
with extended Au pyramids, and transport properties are 
computed after a self-consistent, single-point DFT calcu- 
lation. 



IV. THEORETICAL RESULTS 

In this section we discuss in detail the effect of different 
contact geometries on the molecular conformation of the 
BPDTs and their conduction properties. 
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Figure 5: (Color online) Comparison of the torsion angle ip 
as determined from X-ray measurements and from DFT cal- 
culations in the gas phase (isolated; SR=H in Fig. [l|as well 
as in the molecular junctions (HH, BB, TT; see Fig. |4|). 
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Figure 6: (Color online) Calculated conductance as a function 
of cos'^ for the three types of contact geometries HH, BB, 
and TT. The angles p> are those of the computed junction 
geometries. In all cases, lines show best fits for G = acos^ ip 
with flHH = 2.3-10"^Go (solid), obb = 1.2-10"^Go (dashed), 
aTT = 1-1 • lO'^Go (dash-dotted). 



Molecular Conformation 



B. Conductance 



Fig. [5] shows the torsion angle between the two phenyl 
rings for the molecules as determined by X-ray measure- 
ments and by DFT calculations in the gas-phase as well 
as in the junction geometries. We notice that gas-phase 
angles (with SR=H in Fig. [ij generally coincide well 
with the angles from the X-ray measurements.^^ The 
discrepancy for M6 by roughly 10° has been observed 
previously.^ It is likely due to differences between gas- 
phase and crystal structures caused by the limited stabi- 
lization of the conformation, when there is just a single 
methyl group on each phenyl ring. Note that no X-ray 
structure measurement exists for MO. 

For the contacted molecules deviations of (p from the 
gas-phase conformation are small for HH, but can be 
larger for the BB and TT geometries. This is expected 
from the discussion in Sec. lHI Cl The conformation of the 
alkyl-bridged BPDTs Ml to M5 is very stable. A slight 
trend of increasing tp variations for the molecules with 
the longer, configurationally more flexible alkyl chains 
can be recognized, however. The torsion angles of MO 
and M6 result from the balance between conjugation and 
modest steric repulsion effects due to H atoms or single 
CH3 groups in the ortho position with respect to the ring- 
connecting carbons. Therefore, their (p should be rather 
sensitive to the geometric constraints in the contacts or 
the charge transfer between the molecule and the elec- 
trodes. As a result, deflections of ip from the gas phase 
values of up to 40° occur in the calculations. In con- 
trast, the additional methyl side groups in M7 efficiently 
stabilize ip^ 



In Fig. [6] we present the computed conductance val- 
ues as a function of the torsion angle (p, which the 
biphenyl molecules adopt in the optimized junction ge- 
ometries. On the linear conductance scale we find a rea- 
sonable G — a cos^ (p dependence for all binding situ- 
ations with best fit coefficient^^ ohh = 2.3 • 10~^G'o, 
ABB = 1-2 • lO^^Go, axT = 1-1 • lO^^Go. This behavior 
is characteristic for off-resonant charge transport dom- 
inated by TT-TT coupling and is consistent with the ex- 
perimental observations. Fig. [7] shows, for the sample 
molecule M2, that irrespective of the coordination site, 
the transport is indeed off-resonant and dominated by 
the highest occupied molecular orbital (HOMO) level. 

The results suggest large variations in the conduction 
properties for the different coordination sites of sulfur to 
gold. The conductance of junctions with HH geometry is 
roughly one order of magnitude lower as compared to BB 
and TT, with the sequence of slopes ohh ^ att ~ asB- 
Similar behavior of the conductance of dithiolated aro- 
matic molecules on binding site has been reported before 
by other authors. ^"^ 

For aliphatic alkane molecules, the conductance in the 
bridge-bonded configuration was reported to be higher 
than in the top-bonded one.^^ While these findings are 
compatible with our res ults in Fig. [6] transport through 
alkanes is cr-lik^i^ElIiZl and hence differs substantially 
from the typical 7r-dominated transport through aro- 
matic molecules. Beside the coordination of the an- 
choring group the molecular tilt, which determines the 
overlap of the delocalized tt electrons with the electrode, 
hence plays a crucial role for the conductance of aromatic 
molecules. ''^ We discuss these aspects further below. 
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Figure 7: (Color online) (a) Sketch of the 7r-orbital TBM 
used to describe the transport through BPDT molecules, eo 
is the onsite energy, identical for all carbon atoms, t the cou- 
pling between nearest-neighbor atoms on each phenyl ring, 
t' = t cos if the inter-ring coupling, and ip the torsion angle 
realized in the particular junction geometry. The terminal 
atoms of the biphenyl backbone to the left and right are in- 
dexed a and uu, respectively, and their couplings to the L, R 
electrodes are described by the self energies E£,E5^. (b-d) 
Transmission as a function of energy for M2 in the different 
junction geometries HH, BB, and TT. The solid line is the 
DFT result, the dashed line the fit with the TBM, and the 
dash-dotted line corresponds to the LM. Vertical dashed lines 
indicate, in the order of increasing energy, the position of the 
HOMO, the Fermi energy, and the LUMO. 



However, we note that our contact geometries do not al- 
low us to clearly separate the effects of coordination site 
and tilt, since both are changed simultaneously. 

With regard to absolute values, we observe that the 
calculated conductances are three (BB and TT geome- 
tries) and two (HH geometr y) ord ers of magnitude higher 
than the experimental ones.'^l^ We attribute this over- 
estimation mostly to the interpretation of Kohn- Sham 
orbitals as approximate quasi-particle energies .'^^'^ How- 
ever, also the experimentally measured conductances are 
subject to uncertainties. Indeed, we compare our results 
to the "typical" experimental values, as given by the peak 
position in a conductance histogram, and these peaks are 
rather broad. Further variations of molecular conduc- 
tance, for example due to interactions of the molecules 
with the solvent and other close-by biphenyl molecules 
or the influence of vibrations, have not been accounted 
for in our calculations of static junctions in vacuum. The 
differences on a quantitative level remain as a major chal- 
lenge for future work. 

Finally, we note that our calculations do not repro- 
duce the experimental deviation observed for Ml.'^ Since 
If is unchanged upon contacting (see Fig. [s]), we can ex- 
clude an explanation based on conformational changes, 
which would decrease the degree of conjugation and lead 
to reduced conductances. In spite of the slightly bent 
structure due to the short CH2 bridge (see Fig. [T|), the 
intact Ml shows the highest calculated conductances for 
all coordination geometries (see Fig. [6|. 



C. Analysis of transmission resonances 

In order to understand better the charge transport 
through the BPDT single-molecule junctions, we ana- 
lyze the transmission in terms of a TBM and a LM. We 
use the TBM of Ref. [3SJ which is sketched in Fig. ^a,). 
It describes the delocalized 7r-electron system, relevant 
for transport away from the perpendicular orientation 
of the phenyl rings. The Hiickel-like, molecular Hamilto- 
nian contains three parameters, namely the onsite energy 
eo of each carbon atom, the hopping t between nearest- 
neighbor atoms on each ring, and the torsion angle if, 
specific to the considered junction geometry. Together, 
t and ip determine the matrix element between the ring- 
connecting carbon atoms t' = t cos ip. For the description 
of transport we make use of the wide-band approxima- 
tion, according to which the retarded self energy due 
to the coupling to the electrode X = L, R is determined 
by the line-broadening matrix Tx as = — iPx/Z (cf. 
Eq. A6). We assume a symmetric junction F = {TL)aa = 
{^rKju) and, in line with the nearest-neighbor coupling 



in the molecule, consider the self energy to be non van- 
ishing only on the terminal carbon atoms a and uj of the 
biphenyl backbone [see Fig. [Tf^a)]. The TBM is hence 
characterized by the four parameters (/s, eo,f, F, where ip 
is fixed by the considered junction geometry. 

The parameters of the LM are derived from those of 
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eo 
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~AI2 


~M2 


pAf2 


HH 


-4.40 


-2.30 


0.70 


-1.05 


2.25 


0.11 


BB 


-4.02 


-1.95 


1.10 


-0.42 


2.38 


0.18 


TT 


-4.00 


-1.90 


0.96 


-0.36 


2.36 


0.15 



Table I: Parameters of the TBM eo,t, F obtained by fitting 
the DFT-based t{E) curves for M1-M4. The parameters 
Jf^^gf^f"^ of the LM are those derived from the TBM 
for M2. All values are given in units of eV. 



the TBM. For that purpose we solve the non-Hermitian 



eigenvalue problem J2ki^ + = '^M^fc ^^'^ 

lect the complex eigenvalues corresponding to the HOMO 
and the lowest unoccupied molecular orbital (LUMO). In 
the eigenvalue equation Hjk and {^^)jk = (S£ -I- 
are the Hamiltonian matrix and the self-energy matrix of 



the TBM, respectively, and 



■ «7m- 



We measure 



the real part of the complex eigenvalues with respect to 
the Fermi energy, introducing en = ^homo ^ Ep and 
El = eLUMO — Ep- Due to the symmetries of the TBM 
we find for the imaginary parts ^homo = IlumOi £md 
set r = \lHOMo\- From the relation between the TBM 



and the LM discussed in Appendix p3] [see Eq. (B2|], we 
can identify F with the width of the Lorentzian transmis- 
sion resonances related to the HOMO and the LUMO. 
Finally, we determine the transmission for the LM via 
Eq. ( B2 1 as a sum over these two frontier orbitals only. 



The LM is thus characterized by e//, e^, F and is specific 
to a certain molecule and junction geometry, as described 
by ip,eo,t,T in the TBM. 

Using the physically motivated TBM, we have fitted 
the transmission t{E) of the well-conjugated molecules 
M1-M4, as determined by the DFT calculations. Setting 
Lp to the value of the torsion angle realized in the partic- 
ular junction geometry, we place special emphasis on a 
good fit in the region of the HOMO-LUMO gap. When 
such a fit is too ambitious due to the simplicity of the 
TBM, we describe well at least the region between the 
HOMO and Ep, to determine effective parameters for the 
dominant transmission resonance, as well as the position 
of the LUMO peak. In this way, we obtain the values for 
eo,^,r given in Table |T] Specific LM parameters for M2 
are provided in the same table, and DFT, TBM, and LM 
transmission curves for M2 are shown in Fig. [7|^b-d) . 

The differences between the curves of the TBM and the 
LM in Fig.|7];b-d) in the region of the HOMO-LUMO gap 
illustrate approximations related to the neglect of inter- 
ference effects in the LM. Indeed, we find that the trans- 
mission is slightly overestimated when it is regarded as 
the superposition of incoherent transmission resonances. 
In the following we restrict our discussion to the parame- 
ters of the LM for M2, since they are easy to interpret and 
those of the generic TBM contain similar information. 

The data in Table |l] shows very similar values of f 
for the different junction geometries. While the increas- 
ing linewidth r*^^ -^j^en going from TT to BB is con- 
sistent with the expectation of a better electronic cou- 



pling for a higher coordination of the sulfur atom, also 
the molecular tilt plays a role. The perpendicular ori- 
entation of the BPDTs for geometry HH thus leads to a 
reduced T^'^^ . As an important conclusion, the values of 
e^^2 and ef,^^ show that the HOMO is closer to Ep than 
the LUMO by more than 1 eV. In addition, the reduced 
conductance for HH in Fig. [6] is explained by the HOMO 
level being around 0.5 eV further away from Ep than for 
BB and TT. 

We attribute the shift of the HOMO level towards 
lower energies for increasing coordination number of 
the sulfur atoms to the different amounts of trans- 
ferred charge at the molecule-Au interface. Indeed, both 
Lowdin and electrostatic-potential-derived charges yield 
a leakage of electrons from the molecule, including the S 
atoms, to the Au electrodes, when going from TT over 
BB to the HH geometry. Variations of the conductance 
therefore mostly arise from changes in the alignment of 
the HOMO level with respect to the Fermi energy of the 
Au electrodes, and originate from charge redistributions, 
which are sensitive to the coordination site of the sulfur 
atom at the molecule-electrode interface. 



D. Transmission eigenchannels 

To explore further the electron transport through 
BPDT molecules, especially for the situation ~ 90° 
where the TBM looses its validity, we consider the eigen- 
channel decomposition of the conductance and the cor- 
responding wavefunctions. The results are displayed in 
Figs. |8] and [9| 

We observe one dominant eigenchannel, whose trans- 
mission probability is decreasing gradually with increas- 
ing torsion angle for geometries with (p < 80° (Fig. |8|. 
The wavefunction of this channel is formed from those 
p orbitals of the C atoms, which are perpendicular to 
the phenyl-ring planes (see the results for Ml and M2 
in Fig. |9|. The resulting tt orbitals, which comprise the 
terminal thiol groups, hence exhibit nodes in the ring 
planes. The findings agree with the expectation that for 
the small torsion angles, resulting in a high degree of 
conjugation, electric transport should occur via the delo- 
calized 7r-electron system of the BPDTs. 

The TT-TT coupling between th e rings i s suppressed for 
ip « 90°, since it varies as cos (^1121201251 jj^ ^j^jg ^^^se the 
molecular states become more localized on the individual 
rings. The incoming Bloch waves from the leads can still 
couple through the sulfur linker atom into the 7r-electron 
system of one of the rings, but they are back-refiected at 
the ring-connecting carbon atom. This results in a large 
suppression of the transmission (Fig. |8|, and becomes 
manifest in a low amplitude of the wavefunction on the 
second ring (see the results for M5 and M7 in Fig. |9|. 
In this re gime, the tt-ct coupling, proportional to sinip, 
dominatesl^nill! The a character of the wavefunctions is 
apparent from the absence of nodal planes in the phenyl 
ring planes and the high amplitude of the eigenchannel 
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Figure 8: (Color online) Calculated conductance G and the 
conductance Gn = Got„{Ef) with n = 1, ... ,4 of the four 
transmission eigenchannels with the highest contribution to 
G for the set of BPDTs in (a) HH, (b) BB, and (c) TT con- 
figurations. 



wavefunction on the axis which connects the neighboring 
carbon atoms. 

The isolated biphenyl molecules MO and M7 (SR=H in 
Fig. [T]) with if set to 90° possess Z?2d symmetry. Then. 
CT-TT and TT-fT orbitals are degenerate, which should lead to 
two dominant transmission eigenchannels with the same 
contribution to the conductancc^^ftizs, jjQ.^gygj.^ ^^le pres- 
ence of the electrodes generally leads to a low symme- 
try of the junction as a whole and may also modify the 
molecular geometry. Hence, it is interesting to analyze 





HH 


BB 


TT 


M7 


0.95 


0.15 


1.0 • 10"^ 



Table II: Ratio G2/G1 of the highest eigenchannel contribu- 
tions to the conductance for M7 in the three junction geome- 
tries studied. 



the degeneracy of eigenchannels in the different coordi- 
nation geometries for M7 with the nearly perpendicular 
gas phase torsion angle. Using the ratio of the channel 
conductances G2/G1 with G„ = GoTn{Ep) as a measure 
for the degeneracy, we find the values given in Table [iTl 

The data in Table [H] demonstrates the general absence 
of the channel degeneracy and a high sensitivity of G2 / Gi 
to the junction geometry. Only for the HH contact geom- 
etry we find a nearly perfect degeneracy of the two dom- 
inant transmission eigenchamrels. Consistent with this. 
Fig. [9] demonstrates that their wavefunctions are indeed 
of TT-cr and ct-tt type. The degeneracy can be explained by 
the fact that M7 in the HH geometry stands perpendicu- 
lar to the electrodes. The torsion angle of the contacted 
molecule is hence close to those in the gas phase (see Fig. 
|5|, and the overlap of the molecular tt orbitals with the 
electrode states is such that the degeneracy of molecular 
orbitals is not strongly lifted. Therefore, the transmis- 
sion reflects symmetry properties of the molecule. For 
the BB and TT junctions, deviations from the channel 
degeneracy result from the geometric constraints set by 
the electrodes, which cause Lp in the junctions to deviate 
from 90° , and from the asymmetric overlap of the molec- 
ular TT states with the electrode states to the left and 
right. Thus, our results clearly show that the reduced 
symmetry of the complete junction has to be considered 
for transport and not just the symmetry of the isolated 
molecule alone.'^ 

These findings suggest that measurements of the trans- 
mission eigenchannel degeneracy may serve as a seirsi- 
tive probe to determine the coordination geometry in 
biphenyl- type single-molecule junctions. However, there 
are several factors not included in our idealized treat- 
ment. Thus, it would be interesting to study, how 
strongly a finite bias voltage will lift an existing n-a and 
CT-TT channel degeneracy by breaking of the left-right sym- 
metry. Furthermore, also dynamic effects due to vibra- 
tional modes and Jahn- Teller distortions should lead to 
an effective splitting of the two dominant eigenchannels. 
Beside these issues, it remains an experimental challenge 
to determine the coirduction eigenchannel transparencies 
for contacts with a low transmission, since the existing 
techniques, employing superconducting electrode^^^ or 
shot noise yield very low signals in such situations. 

Our calculations illustrate that the alkyl chains do 
not participate significantly in transport, as expected 
from the large gaps between HOMO and LUMO levels of 
alkanes.Ii^ Considering the dominant transmission eigen- 
channels in Fig. [9] we see that there is indeed practically 
no weight of the wavefunction on the alkyl chain, even 
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Figure 9: (Color online) Wavefunction of the dominant, left-incoming transmission eigenchannel for selected BPDT molecules 
in the HH, BB, and TT geometries. The same isosurface value of the wavefunctions is used in all the plots to allow for their 
comparison. However, the isosurface value has been reduced by a factor of 4 on the right phenyl ring of M7 for HH to visualize 
the TT-o" and a-ir character of the two eigenchannel wavefunctions, which yield the same contribution to the conductance. 



for the short chains present in Ml and M2. 



V. CONCLUSIONS 

Motivated by recent experiments,'^!^!] have pre- 
sented a detailed theoretical analysis of the charge trans- 
port properties of Au-BPDT-Au single-molecule junc- 
tions. The three different types of contact geometries 
in our DFT-based study differed in essential aspects at 
the molecule-metal interface. They were mainly the co- 
ordination site of the anchoring sulfur atoms and the tilt 
of the molecule with respect to the electrodes. Given the 
extensive statistical analysis in the experiments, this set 
of geometries is clearly very limited. Furthermore, with- 
out an analysis of the junction formation process, it is 
difficult to make a statement on the probability of their 
occurrence. However, we hope that they can be used to 
describe general trends, such as the influence of molec- 
ular conformation on conductance and the variability of 
transport properties with contact geometry. 

We have investigated electrode-induced changes of the 
molecular conformation due to charge transfer and geo- 
metric constraints and find that they are rather small for 
most molecules and types of junctions considered here. 
Compared to the somewhat larger variations for MO and 
M6, whose (p is not fixed by an alkyl strap or strong steric 
effects, our calculations show that the appropriate design 
of the side groups can help to stabilize the torsion angle. 

The transport calculations confirm a cos^ (p de- 



pendence of the conductance for the well-conjugated 
molecules in each type of junction geometry. This is 
in accordance with the experimental observations and 
is characteristic for off-resonant transport through the 
TT-electron system."^ For biphenyl molecules with torsion 
angles close to the perpendicular orientation, however, we 
observe systematic deviations in our experimental data 
from the cos^ (p law predicted by a simplified 7r-orbital 
TBM. In that regime of a broken conjugation, our anal- 
ysis of DFT-based transmission eigenchannel wavefunc- 
tions reveals residual conductance contributions from a 
pair of 7r-(T-type conduction channels. 

Finally, our calculations suggest that molecular junc- 
tions with sulfur atoms bound to the "hollow" site of gold 
electrodes could exhibit an order of magnitude smaller 
conductance as compared to junctions with sulfur atoms 
bound via "top" or "bridge" sites. Our analysis shows that 
the transport is dominated by the molecular HOMO level 
in all cases, and variations of the conductance arise from 
changes in the alignment of that level with respect to the 
Fermi energy of the Au electrodes. These changes in turn 
originate from the charge transfer between the molecule 
and the electrode, which is sensitive to the coordination 
site of the sulfur atom. 
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Appendix A: Determination of transmission 
eigenchannels 

In this appendix, we provide further details on how 
we determine the transmission eigenchannels, in particu- 
lar their wave functions. The result is equivalent to that 
of Ref. HH However, our procedure avoids the Lowdin 
orthogonalizations and uses, instead, a consistent for- 
mulation in terms of nonorthogonal basis states. This 
reduces the numerical effort and eliminates possible nu- 
merical instabilities resulting from the forward and back- 
ward Lowdin transformations. Isosurfaces of eigenchan- 
nel wavefunctions employing this scheme are plotted in 
Fig.|9] 

To compute the charge transport, we divide the 
nanocontact into a left (L), central (C), and right (i?) 
region (see Fig. |3| and classify the states of our lo- 
cal, nonorthogonal basis je^) accordingly. The C region 
is assumed to be long enough to neglect the elements 
Hlr = Hjjj^ and Slr = S^j^ of the real and symmetric 
Hamiltonian Hjk = {ej\H\ek) and overlap Sjk = {ej\ek)- 
Here H is the Hamiltonian operator in the combined 
L,C,R space. 

Adopting a notation along the lines of Refs. 1431451 we 
express the energy-dependent transmission as 

TiE)=TT[ALiE)TniE)], (Al) 

where we define the spectral function 

AxiE) - G^cciE)TxiE)G1jciE). (A2) 

Here and below, X stands for either L or R {X — L,R). 
Ax{E) is the contribution to the full spectral density of 
C from scattering states originating in lead X.'^^ In the 
expression. 



where 

g'xxiE) = [{E + ill) Sxx " Hxx]'^ (A5) 
is the retarded Green's function of region X. The matrix 



rx{E)^-2Ira[^^x{E)] 



(A6) 



is the line-broadening matrix. Note that both matrices 
Ax{E) and TxiE) are positive-semidefinite. For nota- 
tional convenience, we will henceforth suppress the en- 
ergy dependence of the quantities. 

In the following we use a basis-independent notation 
with operators such as A]^ and Tji, which are defined by 
their matrix elements in the C space. We also assume 



the existence of a dual basis 



satisfying 



and 1 = E,ec|e'')(ejl- Now let {Al} 



jk 



{e^AL\e 



and (Tii)jk = {ej\Tii\ek)- The matrix elements {TR)jk 
are the components of Eq. ( |A6[ ). They are "covari- 
ant", since the factors Hj^ — ESjk in Eq. (A4| are co- 
variant. The elements (^l)-''"' are also just the com- 
ponents of Eq. (A2). However, they are "contravari- 



ant", since the Green's functions G^q and Gqq [see Eq. 



(|A3|] are defined as the inverse of covariant matrices, i.e. 

± ^v)SJk + Me'^IG'-'le') = S,i. 
Motivated by Eq. ( Al ) we also define the transmission 
probability operator 



Ti = AlTr. 



(A7) 



We will now show how the eigenchannel wavefunctions for 
waves coming in from the left are conveniently obtained 
from the right eigenvectors of Ti in the nonorthogonal 
basis by slightly reformulating the procedure presented 
in Ref. 45, 

Consider the eigenvalue equation 



(A8) 



The eigenvectors \xj) of the Hermitian operator Al are 
orthonormal {{xj\Xk) — Sjk)- Using them we define the 

states \xj) = \f^\Xj) ^^nd the corresponding dual ones 

IX^') = V^IXj) (for all A, ^ 0) so that (x^|xfc) = 5,^- 

It was shown in Ref. |45]that the states \xj) are the de- 
vice part of orthogonal linear combinations of energy- 
normalized scattering states, describing waves coming in 
from the left lead. The transmission eigenchannels |(/)„) 
can be expanded as 



I0n) 



(A9) 



G},c{E) - [{E + ^r7) Sec - Hoc - ^1{E) - ^"-RiE)]-' 

(A3) 

is the retarded Green's function of the device (or G re- 
gion), with 77 > an infinitesimal constant, and G^^ = 
(Gpp)^ is the advanced function. For we need the 
self energies 



j:^xiE) - (Hex - EScx) g'xxiE) {H 



xc 



ESxc) , 
(A4) 



with Cjn = (F I 

may be found from the eigenvalue equation 



The coefficients in this expansion 

(AlO) 



Multiplying by {x^ \ on the left and using Eqs. ( A8 1-( A9 1 
this results iri^ 



k 



j|i R\Xk)Ckn = TnCjn- 



(All) 



11 



Employing the normalization condition X^feCfc^c/ci = Sji, 
the linear combination of Eq. (A9) preserves the en- 
ergy normalization of the left-incoming states. For 
(XjiffllXm) = J2k.i^*kj(^R)kidhn the coefficients dim = 

are determined from Al by multi- 




plying Eq. (A8l by (e^ | on the left and inserting 1. This 



leads to the generalized eigenvalue problem 



Slmdmn 



(A12) 



k.l' 



with dmn = (e"|x„) and J2k.i'^*k]'^kidim. = Sjm- 

Putting these results together, the explicit form of the 
eigenchannel wavefunction for region C in terms of the 
basis functions is obtained from 



) djkCkji 



(A13) 



Using Eqs. ([A3|-(|A6| and ( |Alip -( |AT3| the transmission 
eigenchannel wavefunction can be computed without re- 
sorting to a Lowdin transformation. 



We note that the eigenvalues r„ of Ti in Eq. ( A 10 ) are 
real, since Tn is a Hermitian operator [see Eq. (Alll]. 



It is also easy to show that they agree with the eigen- 
values of more symmetric, Hermitian transmission op- 
erat ors of the form T2 = iP . Given and t„ from 

Eq. ( AlO I and assuming, e.g., t = ^/Tj^G^ y/T^ , the states 

I'P'n) — vT/jl^n) are eigenstates of T2 with the eigenval- 
ues — Tn- Furthermore, it is easy to provd^that the 
eigenvalues satisfy < ?■„ < 1, as expected for transmis- 
sion probabilities. 



Appendix B: Relation between the tight-binding and 
Lorentz model 



The Lorentz model (LM) is frequently used to describe 
the transmission in the field of molecular electronics. 
Typically, a Lorentzian function is fitted to the resonance 
dominating the transmission at the Fermi energy. Here, 
we discuss, in which situation the LM coincides with the 
TBM. 

We consider the non-Hermitian eigenvalue problem 
(H -I- Y7y\ii) = A^l/i) with = + ^7^, the symmet- 
ric and Hermitian Hamilton operator iJ, and the sym- 
metric, but non-Hermitian retarded self-energy operator 



Y7 = S£ + composed of contributions from the L 
and R electrodes. By {Ji\ we denote the left eigenstate 
with the same eigenvalue as the corresponding right 
eigenstate |/^), i.e. {fi\{H + S'') = A^(/2|. The C re- 
gion is assumed to be identical to the molecule in the 
TBM [see Fig. |7ja)]. Using the spectral decomposition 
of the Green's function in the expression for the energy- 
dependent transmission t{E) [see Eqs. (All and (A2|], 
we obtain 



{fl\TLW){i>\tB.\fi) 



[E- 



(Bl) 



where the sum is over all those eigenstates \fi) of the 
biphenyl which obtain a finite linewidth 7^ 7^ by the 
coupling to the electrodes and which hence contribute to 
the transport. 

Let us now make the wide-band approximation 
and consider the energy-independent expression S'' = 
— i(rL+rfl)/2 to be a small perturbation. Within lowest- 
order perturbation theory we obtain A^ = ej^ + 47^ with 

= e°|Ai") and 7^ = (^0|S''|^0). Additionally, we 
assume a symmetric coupling (ri)^^ — {Tfi)i^i^ = F, 
where the indices a, lu refer to those atoms of the biphenyl 
backbone which are closest to the L, R electrodes [see 
Fig. ^B,)] and where local basis states are understood to 
be orthogonal in the spirit of the Hiickel approximation. 
By exploiting the inversion symmetry of the TBM, it fol- 
lows that 7^ = (/i°|fx|At°) since = i, MHM = H, 
and MTlM = Tji with the operator M — describing 
the inversion of the molecule. 

The perturbation theory is valid in the regime T <^ t, 
where t determines the separation between the resonance 
energies relevant for transport. When they are well 
separated, the largest contributions to the transmission 
in Eq. (Bl I arise when ^ — v, since cross-terms are sup- 



pressed by a large off-resonant denominator. In this case 
the transmission is well represented as the sum of inco- 
herent Lorentz resonances 



r{E) « E 



[E 



(B2) 



and the TBM simplifies to the LM. 
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